Spitting in the wind?-The challenges of RNA sequencing for biomarker discovery from saliva.

Forensic trace contextualization, i.e., assessing information beyond who deposited a biological stain, has become an issue of great and steadily growing importance in forensic genetic casework and research. The human transcriptome encodes a wide variety of information and thus has received increasing interest for the identification of biomarkers for different aspects of forensic trace contextualization over the past years. Massively parallel sequencing of reverse-transcribed RNA ("RNA sequencing") has emerged as the gold standard technology to characterize the transcriptome in its entirety and identify RNA markers showing significant expression differences not only between different forensically relevant body fluids but also within a single body fluid between forensically relevant conditions of interest. Here, we analyze the quality and composition of four RNA sequencing datasets (whole transcriptome as well as miRNA sequencing) from two different research projects (the RNAgE project and the TrACES project), aiming at identifying contextualizing forensic biomarker from the forensically relevant body fluid saliva. We describe and characterize challenges of RNA sequencing of saliva samples arising from the presence of oral bacteria, the heterogeneity of sample composition, and the confounding factor of degradation. Based on these observations, we formulate recommendations that might help to improve RNA biomarker discovery from the challenging but forensically relevant body fluid saliva.

Selecting mRNA markers in blood for age estimation of the donor of a biological stain.

RNA has gained a substantial amount of attention within the forensic field over the last decade. There is evidence that RNAs are differentially expressed with biological age. Since RNA can be co-extracted with DNA from the same piece of evidence, RNA-based analysis appears as a promising molecular alternative for predicting the biological age and hence inferring the chronological age of a person. Using RNA-Seq data we searched for markers in blood potentially associated with age. We used our own RNA-Seq data from dried blood stains as well as publicly available RNA-Seq data from whole blood, and compared two different approaches to select candidate markers. The first approach focused on individual gene analysis with DESeq2 to select the genes most correlated with age, while the second approach employed lasso regression to select a set of genes for optimal prediction of age. We present two lists with 270 candidate markers, one for each approach.

Development and validation of an mRNA-based multiplex body fluid identification workflow and a rectal mucosa marker pilot study.

Molecular identification of body fluids and tissues is crucial in order to understand the circumstances of crimes. For that reason, molecular investigations used to identify body fluids/tissues have increasingly been examined recently. Various studies have proved that messenger RNA (mRNA) profiling is a sensitive and robust method for body fluid/tissue identification. The forensically relevant body fluids/tissues blood, semen, saliva, vaginal secretion, menstrual blood and skin have all been detected successfully by applying suitable mRNA assay. However, rectal mucosa, which can be found as evidence in sexual assault cases, has been neglected in forensic investigations. So far there is no mRNA marker to detect rectal mucosa, although anal penetration occurs in a large number of sexual assaults (23.2% of female victims and 50% of male victims). In this study, specific and sensitive mRNA markers for forensically relevant body fluids were adapted and validated in an mRNA multiplex assay for routine casework. This included the implementation of a DNA/RNA re-extraction method for automated extraction that can be integrated into casework without loss of DNA. This re-extraction method and the mRNA multiplex assay were tested using casework samples. PCR-primers were designed for the identification of rectal mucosa and the more effective marker MUC12 was integrated into an extended multiplex assay. The result of our study is a highly specific and sensitive mRNA multiplex assay plus an automated DNA/RNA re-extraction method, that can be integrated into casework and identify rectal mucosa for the first time.

mRNA profiling of mock casework samples: Results of a FoRNAP collaborative exercise.

In recent years, forensic mRNA profiling has increasingly been used to identify the origin of human body fluids. By now, several laboratories have implemented mRNA profiling and also use it in criminal casework. In 2018 the FoRNAP (Forensic RNA Profiling) group was established among a number of these laboratories with the aim of sharing experiences, discussing optimization potential, identifying challenges and suggesting solutions with regards to mRNA profiling and casework. To compare mRNA profiling methods and results a collaborative exercise was organized within the FoRNAP group. Seven laboratories from four countries received 16 stains, comprising six pure body fluid / tissue stains and ten mock casework samples. The laboratories were asked to analyze the provided stains with their in-house method (PCR/CE or MPS) and markers of choice. Five laboratories used a DNA/RNA co-extraction strategy. Overall, up to 11 mRNA markers per body fluid were analyzed. We found that mRNA profiling using different extraction and analysis methods as well as different multiplexes can be applied to casework-like samples. In general, high input samples were typed with high accuracy by all laboratories, regardless of the method used. Irrespective of the analysis strategy, samples of low input or mixed stains were more challenging to analyze and interpret since, alike to DNA profiling, a higher number of markers dropped out and/or additional unexpected markers not consistent with the cell type in question were detected. It could be shown that a plethora of different but valid analysis and interpretation strategies exist and are successfully applied in the Forensic Genetics community. Nevertheless, efforts aiming at optimizing and harmonizing interpretation approaches in order to achieve a higher consistency between laboratories might be desirable in the future. The simultaneous extraction of DNA alongside RNA showed to be an effective approach to identify not only the body fluid present but also to identify the donor(s) of the stain. This allows investigators to gain valuable information about the origin of crime scene samples and the course of events in a crime case.

Inter-laboratory study on standardized MPS libraries: evaluation of performance, concordance, and sensitivity using mixtures and degraded DNA.

We present results from an inter-laboratory massively parallel sequencing (MPS) study in the framework of the SeqForSTRs project to evaluate forensically relevant parameters, such as performance, concordance, and sensitivity, using a standardized sequencing library including reference material, mixtures, and ancient DNA samples. The standardized library was prepared using the ForenSeq DNA Signature Prep Kit (primer mix A). The library was shared between eight European laboratories located in Austria, France, Germany, The Netherlands, and Sweden to perform MPS on their particular MiSeq FGx sequencers. Despite variation in performance between sequencing runs, all laboratories obtained quality metrics that fell within the manufacturer's recommended ranges. Furthermore, differences in locus coverage did not inevitably adversely affect heterozygous balance. Inter-laboratory concordance showed 100% concordant genotypes for the included autosomal and Y-STRs, and still, X-STR concordance exceeded 83%. The exclusive reasons for X-STR discordances were drop-outs at DXS10103. Sensitivity experiments demonstrated that correct allele calling varied between sequencing instruments in particular for lower DNA amounts (≤ 125 pg). The analysis of compromised DNA samples showed the drop-out of one sample (FA10013B01A) while for the remaining three degraded DNA samples MPS was able to successfully type ≥ 87% of all aSTRs, ≥ 78% of all Y-STRs, ≥ 68% of all X-STRs, and ≥ 92% of all iSNPs demonstrating that MPS is a promising tool for human identity testing, which in return, has to undergo rigorous in-house validation before it can be implemented into forensic routine casework.

"The acid test"-validation of the ParaDNA® Body Fluid ID Test for routine forensic casework.

The identification of the cellular origin and composition of crime scene-related traces can provide crucial insight into a crime scene reconstruction. In the last decade, especially mRNA-based body fluid and tissue identification (BFI) has been vigorously examined. Besides capillary electrophoretic (CE) and real-time quantitative PCR (RT-qPCR)-based approaches for mRNA detection, melt curve analysis bears potential as a simple-to-use method for BFI. The ParaDNA® Body Fluid ID Test relies on HyBeacon® probes and was developed as a rapid test for mRNA-based BFI of six different body fluids: vaginal fluid, seminal fluid, sperm cells, saliva, menstrual, and peripheral blood. The herein presented work was performed as an "acid test" of the system and should clarify whether the approach matches the requirements of forensic routine casework in German police departments. Tested samples consisted of single source as well as of mixed samples.

As solid as a rock-comparison of CE- and MPS-based analyses of the petrosal bone as a source of DNA for forensic identification of challenging cranial bones.

Short tandem repeat (STR) typing from skeletal remains can be a difficult task. Dependent on the environmental conditions of the provenance of the bones, DNA can be degraded and STR typing inhibited. Generally, dense and compact bones are known to preserve DNA better. Several studies already proved that femora and teeth have high DNA typing success rates. Unfortunately, these elements are not present in all cases involving skeletal remains. Processing partial or singular skeletal elements, it is favorable to select bone areas where DNA preservation is comparably higher. Especially, cranial bones are often accidentally discovered during criminal investigations. The cranial bone is composed of multiple parts. In this examination, we evaluated the potential of the petrous bone for human identification of skeletal remains in forensic case work. Material from different sections of eight unknown cranial bones and-where available-additionally other skeletal elements, collected at the DNA department of the Institute of Legal Medicine in Ulm, Germany, from 2010 to 2017, were processed with an optimized DNA extraction and STR typing strategy. The results highlight that STR typing from the petrous bones leads to reportable profiles in all individuals, even in cases where the analysis of the parietal bone failed. Moreover, the comparison of capillary electrophorese (CE) typing to massively parallel sequencing (MPS) analysis shows that MPS has the potential to analyze degraded human remains and is even capable to provide additional information about phenotype and ancestry of unknown individuals.

Conserved co-regulation and promoter sharing of hoxb3a and hoxb4a in zebrafish.

The expression of zebrafish hoxb3a and hoxb4a has been found to be mediated through five transcripts, hoxb3a transcripts I-III and hoxb4a transcripts I-II, driven by four promoters. A "master" promoter, located about 2 kb downstream of hoxb5a, controls transcription of a pre-mRNA comprising exon sequences of both genes. This unique gene structure is proposed to provide a novel mechanism to ensure overlapping, tissue-specific expression of both genes in the posterior hindbrain and spinal cord. Transgenic approaches were used to analyze the functions of zebrafish hoxb3a/hoxb4a promoters and enhancer sequences containing regions of homology that were previously identified by comparative genomics. Two neural enhancers were shown to establish specific anterior expression borders within the hindbrain and mediate expression in defined neuronal populations derived from hindbrain rhombomeres (r) 5 to 8, suggesting a late role of the genes in neuronal cell lineage specification. Species comparison showed that the zebrafish hoxb3a r5 and r6 enhancer corresponded to a sequence within the mouse HoxA cluster controlling activity of Hoxa3 in r5 and r6, whereas a homologous region within the HoxB cluster activated Hoxb3 expression but limited to r5. We conclude that the similarity of hoxb3a/Hoxa3 regulatory mechanisms reflect the shared descent of both genes from a single ancestral paralog group 3 gene.

The Trichoplax PaxB gene: a putative Proto-PaxA/B/C gene predating the origin of nerve and sensory cells.

Pax genes play key regulatory roles in embryonic and sensory organ development in metazoans but their evolution and ancestral functions remain widely unresolved. We have isolated a Pax gene from Placozoa, beside Porifera the only metazoan phylum that completely lacks nerve and sensory cells or organs. These simplest known metazoans also lack any kind of symmetry, organs, extracellular matrix, basal lamina, muscle cells, and main body axis. The isolated Pax gene from Trichoplax adhaerens harbors a paired domain, an octapeptide, and a full-length homeodomain. It displays structural features not only of PaxB and Pax2/5/8-like genes but also of PaxC and Pax6 genes. Conserved splice sites between Placozoa, Cnidaria, and triploblasts, mark the ancient origin of intron structures. Phylogenetic analyses demonstrate that the Trichoplax PaxB gene, TriPaxB, is basal not only to all other known PaxB genes but also to PaxA and PaxC genes and their relatives in triploblasts (namely Pax2/5/8, Pax4/6, and Poxneuro). TriPaxB is expressed in distinct cell patches near the outer edge of the animal body, where undifferentiated and possibly multipotent cells are found. This expression pattern indicates a developmental role in cell-type specification and/or differentiation, probably in specifying-determining fiber cells, which are regarded as proto-neural/muscle cells in Trichoplax. While PaxB, Pax2/5/8, and Pax6 genes have been linked to nerve cell and sensory system/organ development in virtually all animals investigated so far, our study suggests that Pax genes predate the origin of nerve and sensory cells.

Comparative genomic analysis of vertebrate Hox3 and Hox4 genes.

We used a comparative genomic approach to identify putative cis-acting regulatory sequences of the zebrafish hoxb3a and hoxb4a genes. We aligned genomic sequences spanning the clustered Hoxb1 to Hoxb5 genes from pufferfish, mice, and humans with the zebrafish hoxba and hoxbb cluster sequences. We identified multiple blocks of conserved sequences in non-coding regions within and surrounding the Hoxb3/b4 gene locus; a subset of these blocks are conserved in the zebrafish hoxbb cluster, despite loss of hoxb3/b4 genes. Overall, we find that the architecture of the Hoxb3/b4 loci and of the conserved sequence elements is very similar in teleosts and mammals. Our analyses also revealed two alternative transcripts of the zebrafish hoxb3a gene and an exon sequence unusually located 10 kb upstream of adjacent hoxb4a; an equivalent murine Hoxb3 exon has not yet been confirmed. We show that many of the Hoxb3/b4 conserved non-coding sequences correlate with functional neural enhancers previously described in the mouse. Further, within the conserved non-coding sequences we have identified binding sites for transcription factors, including Kreisler/Valentino, Krox20, Hox, and Pbx, some of which had not been previously described for the mouse. Finally, we demonstrate that the regulatory sequences of zebrafish hoxa3a are divergent with respect to the mouse ortholog Hoxa3, or the paralog hoxb3a. Despite limited conservation of regulatory sequences, zebrafish hoxa3a and hoxb3a genes share very similar expression profiles.